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1.  Introduction 


With  clever  integration  of  microelectromechanical  sensors  and  actuators  into  new  projectile 
configurations,  future  weapon  systems  hope  to  achieve  a  notable  leap  in  accuracy  and  lethality. 
Because  some  sensors  and  control  mechanisms  require  different  levels  of  spin  to  operate 
properly,  new  multicomponent  projectile  configurations  are  being  investigated.  The  dual-spin 
projectile  configuration  is  one  such  concept  that  consists  of  forward  and  aft  components, 
connected  through  a  bearing,  which  roll  at  different  rates. 

A  potential  control  mechanism  for  smart  weapons  trajectory  control  is  a  cluster  of  lateral 
pulsejets.  Each  pulsejet  imparts  a  relatively  short  duration,  yet  large  lateral  body  force  on  the 
projectile.  By  mounting  a  set  of  pulsejets  on  either  the  forward  or  aft  section  of  a  dual-spin 
projectile,  a  specific  trajectory  can  be  tracked.  When  designing  a  projectile  pulsejet  flight 
control  system  to  minimize  tenninal  miss  distance,  pulsejet  firing  logic  relies  on  real-time 
estimation  of  the  effect  of  firing  a  particular  jet  on  the  miss  distance  at  an  arbitrary  time  and 
system  state.  In  theory,  this  can  be  accomplished  by  numerically  integrating  the  nonlinear 
equations  of  motion  from  a  given  point  to  the  target.  However,  from  a  practical  point  of  view, 
real-time  numerical  integration  of  the  nonlinear  equations  of  motion  from  a  given  state  to  the 
target  at  each  flight  control  system  computation  cycle  represents  too  large  a  computational 
burden  on  the  on-board  microprocessor.  Hence,  computationally-efficient  trajectory  prediction 
of  smart  weapon  projectile  configurations  is  needed  for  insertion  into  control  laws.  Along  these 
lines,  Guidos  and  Cooper  ( 1 )  developed  closed  form  expressions  for  swerve  of  a  rigid  nonrolling 
projectile  with  a  lateral  pulsejet  applied.  Their  analysis  allows  for  finite  duration  impulse  loads 
and  showed  that  under  mild  conditions  on  spin  rate,  singular  and  finite  duration  impulses  yield 
the  same  effect  on  the  trajectory. 

Although  the  dynamics  of  dual-spin  spacecraft  have  been  extensively  studied  (2-12),  work 
on  the  atmospheric  flight  mechanics  of  dual-spin  projectiles  has  been  comparatively  sparse. 
Smith  et  al.  (13)  developed  a  dynamic  model  of  a  dual-spin  projectile  and  employed  this  model 
to  simulate  a  smart  artillery  shell.  More  recently,  Costello  and  Peterson  (14)  reported  a  more 
general  dual-spin  projectile  dynamic  model  which  permits  asymmetric  mass  properties  and 
generalized  bearing  friction.  This  model  was  subsequently  used  to  develop  a  linear  theory  and 
stability  boundaries  for  this  configuration.  The  work  reported  here  extends  the  dual-spin 
projectile  linear  equations  to  account  for  the  effect  of  lateral  pulsejets  in  the  prediction  of  angular 
and  swerving  motion.  The  swerving  dynamics  are  solved  in  closed  form,  resulting  in 
computationally-simple  algebraic  expressions  for  the  projectile  trajectory  under  the  influence  of 
lateral  pulsejets. 
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2.  Dual-Spin  Projectile  Dynamic  Model 


The  dynamic  model  for  a  dual-spin  projectile  admits  a  total  of  seven  degrees-of-freedom, 
including  the  components  of  the  position  vector  of  the  mass  center  and  the  Euler  yaw,  pitch,  aft 
body  roll,  and  forward  body  roll  angles.  Figures  1  and  2  provide  indicate  coordinates  used  to 
specify  position  and  orientation  for  the  dual-spin  projectile  dynamic  model.  Linear  projectile 
theory  is  used  to  simplify  the  highly  nonlinear  dynamic  equations  to  a  set  of  linear  ordinary 
differential  equations.  More  details  on  projectile  linear  theory  can  be  found  in  McCoy  (15).  The 
linear  equations  developed  by  Costello  and  Peterson  (14)  neglected  gravity,  since  their  primary 
interest  was  dynamic  stability,  and  also  did  not  include  lateral  pulsejet  forces  in  the  formulation. 
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Figure  2.  Projectile  orientation  definitions. 


Equations  1-10  remove  these  two  limitations  and  represent  a  modified  set  of  dynamic  equations 
describing  the  motion  of  a  dual-spin  projectile  in  atmospheric  flight. 

V'  =  -  pSDC  xo  v  •  (1) 

2  m 


,,  _D 
YF  ~  y  P F  * 


li 

-  it,  -  ^ 

H 

T 

ill 

Pf 

J Pa . 

1* 

\M\ 

W' 

V' 

-A 

0 

0 

-D 

V 

YjD/(mV ) 

w 

0 

-A 

D 

0 

w 

Dg  IV  +  ZjD  l(m  V) 

t 

q 

BID 

C/D 

E 

-F 

s 

q 

>  r  < 

-R^D/iVIyy) 

r' 

—C  /  D 

BID 

F 

E 

r 

RJXY,D/(V/W) 

(2) 

(3) 

(4) 


(5) 


e ' 


(6) 


¥ 


D 

—  r  • 
V 


(V) 


3 


x'  =  D. 

y  ■  =  —  v  +  «/£> 
F 


z '  =  —  w  -  9D  ■ 
V 


B  = 


pSD 
2  m 


mD 


A  = 


pSD 
2  m 


(CNA)~ 


y  |  + 


C  = 


r 

A, 


-pf+(R 


+  r 

ax  ax 


rA 

^  NP  A 


E  = 


"  pSD~ 
2  m 

mD 

IT 

L  YY  J 

r  • 

^  MA 

pSD  1 

m  D 2 

r 

MQ 

2  m 

IT 

1  YY 

2 

_  D  (I  XX  P  F  +  1  XxP  a  ) 


V 


H  = 


pSD2CFLP 


DC,, 


47  VI 

^ 1  XX  r  1 XX 


DC , 


J= 


V 


VI 


XX 


dc ; 


r  =  - 


F 


7  = 


VI 

XX 

pSD^C^p  DCV 


41 


XX 


VI 


XX 


M  = 


N  = 


pSD2VC^D 

2/fv 


pSD2VC^D 


21 


A 


CMA  =  (Rfx  +  rfx)cFNA+{Rax  +  rax)C; 


IT  =  IF  +m  r  2  +  IA  +  m  r  2 . 

1 YY  1YY  Hffx  1YY  ,Ua'ax  ’ 


(8) 

(9) 

(10) 


In  equations  1-10,  the  ‘  symbol  denotes  differentiation  with  respect  to  nondimensional  arc 
length,  5  .  The  constants  used  in  the  previous  equations  are  defined  as  follows: 
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The  impulse  force  components  Ty  and  Z7  represent  impulse  forces  that  are  applied  to  the 
projectile  in  the  nonrolling  reference  frame.  These  loads  and  the  associated  moment  arm,  RfX  , 
are  effective  forces  and  moments.  Due  to  aerodynamic  interaction,  the  force  and  moment 
generated  by  the  lateral  pulsejets  is  in  general  different  than  what  would  be  measured  on  a  test 
stand.  Also,  for  this  report,  the  impulse  forces  are  mathematically  treated  as  singular  impulses. 
The  bearing  coupling  rolling  resistance  is  assumed  to  be  of  the  hydrodynamic  type. 


3.  Analytic  Solution 


Pulsejet  forces  and  moments  that  appear  in  equation  5  are  assumed  to  occur  over  a  relatively 
short  duration  of  time  compared  to  the  motion  of  the  projectile;  hence,  their  effect  is  impulsive  in 
nature.  An  impulsive  force  acting  on  a  rigid  body  creates  a  discontinuity  in  the  translational  and 
rotational  velocity  of  the  body  at  the  time  of  the  impulse  (16).  To  simulate  motion  of  the 
projectile  with  pulsejets  acting  on  the  system,  a  piecewise  solution  is  constructed  that  consists  of 
free-flight  and  pulsejet  discontinuity  segments.  Expressions  for  free-flight  solution  segments  are 
given  in  this  section,  while  the  pulsejet  state  discontinuity  conditions  are  provided  in  the  next 
section. 


As  seen  in  equation  1,  the  total  velocity  is  dynamically  uncoupled  from  the  system,  but  appears 
in  a  nonlinear  fashion  in  the  remaining  state  dynamic  equations.  As  is  customary  in  projectile 
linear  theory  (15),  it  is  assumed  that  the  total  velocity  changes  slowly  with  respect  to  the  other 
states  of  the  system  so  it  can  be  treated  as  a  constant  in  the  other  dynamic  equations.  When 
viewed  in  this  light,  the  equations  of  motion  become  linear  and  the  equations  of  motion  for  angle 
of  attack  and  roll  motion  become  uncoupled. 

The  total  velocity  of  the  projectile  exhibits  exponential  decay  as  the  round  flies  down  range. 
Equation  24  provides  the  total  velocity  solution: 

-pSDC  yns/(2m)  .... 

V  =  VQe  xo  ■  (24) 


The  solution  to  the  roll  rate  equations  is  aided  by  first  defining  the  decay  rates,  as  shown  in 
equation  25. 
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With  the  roll  rate  decay  rate  given  by  equation  25,  the  roll  rate  solution  can  be  written  as 
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As  given  in  equation  5,  the  equations  of  motion  for  epicyclic  pitching  and  yawing  behavior 
consists  of  coupled  motion  of  v,  w,  q ,  and  r  .  The  slow  and  fast  modes  of  the  epicyclic 

dynamic  equations  are  given  by  equations  32  and  33: 
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and 


=  -L F-VF2  -4 Cl 


Corresponding  to  the  modes  given  in  equations  32  and  33,  the  right  and  left  Eigenvector  matrices 
are  given  by  equations  38  and  39,  respectively. 
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The  constants  utilized  in  equations  38  and  39  are  defined  as  follows.  Note  that  the  *  symbol 
denotes  complex  conjugate. 

K  =  (E-A)  +  2A  +  iF, 

Q  =  (E-A)2  +4AE  +  4C-F2  +2i(F(E-A)  +  2(AF  +  B)), 

R  =  (E-A)  +  2A-iF, 

S  =  (E  -  A)2  +  4AE  +  4C  -  F2  -2i(F  (E  -  A)  +  2(AF  +  B)) , 

„_kWe) 


7 


The  vector  % ,  which  represents  the  initial  conditions  with  the  steady  state  portion  removed,  is 
used  to  define  the  homogeneous  response  of  the  epicyclic  equations. 
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Using  the  definitions  provided  by  equations  38,  39,  and  46,  the  homogeneous  solution  for  each 
epicyclic  state  variable  is  expressed  as  a  combination  of  the  fast  and  slow  modes. 
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while  the  particular  solution  is  given  by  equation  5 1 . 
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Expressions  for  the  modal  amplitudes  and  phase  angles  used  in  equations  47-50  are  provided  in 
the  appendix.  The  total  epicyclic  solution  is  the  summation  of  the  homogeneous  and  particular 
solutions. 
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The  solutions  for  the  Euler  pitch  and  yaw  angles  are  obtained  by  integrating  the  solutions  for  the 
pitch  rate  and  yaw  rate,  respectively. 
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The  trajectory  or  swerve  solution  is  then  obtained  by  substituting  equations  47,  48,  and  49  into 
equations  9  and  10  and  subsequently  integrating.  The  final  results  are  given  by  equations  57  and  58. 
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Equations  57  and  58 
dual-spin  projectile ; 


<t>SS  +  ewS-Jr)-  sin  {Ps  -x)-</>S  COS  idwS  ~x)s)  ■ 

J 

•  expressions  for  the  altitude  and  cross  range  < 
)th  the  cross  range  and  altitude  solutions  cont 
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constant,  linear,  and  oscillatory  terms.  The  altitude  formula  has  an  additional  parabolic  tenn  to 
account  for  trajectory  bending  due  to  gravity. 


4.  Pulsejet  Force  and  Moment  Conditions 


The  net  effect  of  a  pulsejet  force  is  a  discontinuity  in  the  translational  and  rotational  velocity  of 
the  body  at  the  time  of  the  impulse.  The  discontinuity  conditions  are  found  by  transforming 
equation  5  to  the  Laplace  domain  and  noting  the  additional  impulse  terms  adding  with  the  pre¬ 
impulse  initial  conditions.  The  results  are  given  by  equations  59-62. 


+  ZD 

vo  =  fi>  +  — ' 
mV 

+  Z,D 

™o=wo  + 


mV 


%  =  q0 


IyyV 


(59) 

(60) 

(61) 


=  ro  + 


RgYjD 

IyyV 


(62) 


In  equations  52-55,  the  superscript  “+”  denotes  after  impulse  and  the  superscript  denotes 
before  impulse.  Since  the  epicyclic  solution  amplitude  and  phase  angles  (QvF,  QvS  ,  QvF,  flwF , 

QwS ,  ,  q?s  ,  q,.f ,  q„5  ,  evF ,  evS ,  evF ,  ewF ,  ewS ,  oqF ,  eqS ,  erF  ,  and  0rS  )  all  depend  on  v0 , 

vv0 ,  q0,  and  r0-% ,  the  dynamic  model  must  be  updated  at  s 1  to  accurately  account  for  the 
discontinuity  in  the  epicyclic  solution  fonnulas. 


5.  Results 


In  order  to  investigate  the  utility  of  the  simple  trajectory  fonnulas  derived  previously,  the  results 
that  follow  compare  the  trajectory  of  a  dual-spin  projectile  computed  by  numerically  integrating 
the  fully  nonlinear  equations  of  motion  (13)  with  the  linear  theory  analytical  solution.  The 
numerical  solution  was  generated  using  a  fixed  step  fourth-order  Runge  Kutta  algorithm  with  a 
time  step  of  0.00001  s.  Tables  1  and  2  list  the  physical  parameters  and  state  conditions  for  the 
projectile  considered  in  trajectory  comparison.  The  linear  model  is  updated  for  each  change  in 
range  of  one  caliber.  In  all  cases,  the  numerical  nonlinear  simulation  and  the  analytic  closed 
form  solution  are  in  agreement.  Figures  3  and  4  show  the  swerve  of  the  projectile.  A  difference 
of  approximately  10  ft  between  the  numerical  nonlinear  solution  and  the  analytic  solution  at  a 
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Table  1.  Projectile  physical  properties. 


Parameter 

Value 

Projectile  length 

4.5  ft 

Projectile  reference  diameter 

0.23  ft 

Projectile  weight 

15.9  lb 

Composite  body  mass  center  location 

2.81  ft  (from  base) 

Forward  body  mass  center  location 

3.85  ft  (from  base) 

Aft  body  mass  center  location 

1.35  ft  (from  base) 

Forward  body  roll  inertia 

0.0022  slug  ft2 

Forward  body  pitch  inertia 

0.019  slug  ft2 

Aft  body  roll  inertia 

0.002  slug  ft2 

Aft  body  pitch  inertia 

0.42  slug  ft2 

Bearing  location 

3.18  ft  (from  base) 

Bearing  viscous  damping  coefficient 

0.0001 

Lateral  pulsejet  stationline  location 

3.8  ft  (from  base) 

Lateral  pulsejet  buttline  location 

0.12  ft 

Lateral  pulsejet  waterline  location 

0.0  ft 

Table  2.  Projectile  state  conditions. 


State  Variable 

Initial 

Before  Pulse 

After  Pulse 

X 

0.0  ft 

419.0  ft 

419.0  ft 

V 

-0.76  ft 

-1.0  ft 

-1.0  ft 

z 

-287.0  ft 

-375.0  ft 

-375.0  ft 

(j>F 

4.6° 

O 

OO 

O 

00 

'n 

<t>A 

168.1° 

156.5° 

156.5° 

e 

0.2° 

0.2° 

0.2° 

Y 

-0.0007° 

0.0006° 

0.0006° 

V 

2497.0  ft/s 

2420.0  ft/s 

2420.0  ft/s 

V 

0.17  ft/s 

-0.1  ft/s 

1.8  ft/s 

w 

-0.64  ft/s 

0.4  ft/s 

-0.9  ft/s 

Pf 

7.3  rad/s 

6.7  rad/s 

6.7  rad/s 

Pa 

-10.1  rad/s 

-83.8  rad/s 

-83.8  rad/s 

q 

-0.006  rad/s 

-0.01  rad/s 

2.14  rad/s 

r 

0.001  rad/s 

0.0004  rad/s 

3.04  rad/s 

11 


Figure  3.  Cross  range  vs.  range. 


Figure  4.  Altitude  vs.  range. 
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range  of  3750  ft  is  apparent  in  the  altitude  trace.  When  viewed  from  the  rear,  the  projectile  drifts 
slightly  to  the  left  prior  to  the  pulsejet  firing.  After  the  pulsejet  is  fired,  the  pulsejet  turns  and 
drifts  to  the  right.  Toward  the  end  of  the  trajectory,  the  lateral  pulsejet  has  changed  the  cross 
range  impact  point  by  approximately  17  ft.  Figures  5  and  6  show  that  the  action  of  the  pulsejet 
induces  Euler  angle  changes  exceeding  5°.  Because  of  the  roll  orientation  when  the  pulsejet  was 
fired,  the  altitude  response  does  not  show  an  effect  from  the  pulsejet.  The  fins  mounted  on  the 
aft  body  along  with  inherent  roll  damping  combine  to  create  the  aft  body  roll  rate  history  shown 
in  figure  7.  The  forward  body  roll  rate  history  is  driven  by  the  viscous  friction  in  the  bearing  and 
the  roll  rate  of  the  aft  body.  Figure  8  displays  the  difference  in  the  predicted  impact  point  of  the 
numerical  nonlinear  simulation  and  the  analytic  solution  as  a  function  of  the  linear  model  update 
interval.  At  the  range  of  3750  ft,  the  impact  error  prediction  is  below  10  ft.  Updating  the  linear 
model  used  for  the  analytical  solution  beyond  4000  cal  of  travel  results  in  significant  increase  in 
impact  point  error.  Figure  9  examines  the  potential  change  in  impact  point  location  as  a  function 
of  down  range  position.  After  the  epicyclic  dynamics  have  settled,  the  swerve  solution  is 
dominated  by  the  linear  tenn.  Firing  a  pulsejet  changes  the  slope  of  the  linear  term.  As  the 
projectile  closes  on  the  target,  the  potential  change  in  position  of  the  round  at  the  target  due  to 
firing  a  lateral  pulsejet  decreases,  since  a  reduced  distance  to  the  target  reduces  the  linear  swerve 
changes  at  the  target.  Defining  control  authority  as  the  potential  change  in  swerve  at  the  target, 
the  control  authority  of  a  lateral  pulsejet  decreases  as  the  projectile  travels  down  range. 
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Figure  6.  Yaw  angle  vs.  range. 
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Figure  7.  Roll  rate  vs.  range. 
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Figure  9.  Induced  pulsejet  swerve  vs.  range  to  target. 
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6.  Conclusions 


A  simple  set  of  formulas  has  been  developed  for  the  swerve  motion  of  a  dual-spin  projectile 
under  the  action  of  lateral  pulsejets.  The  solutions  provide  a  computationally  low  cost  and 
reasonably  accurate  means  of  predicting  impact  points  of  a  projectile  that  can  be  employed  in 
future  smart  weapon  flight  control  systems.  For  the  best  impact  point  prediction  performance, 
the  linear  model  that  underlies  the  analytic  solution  must  be  periodically  updated.  For  the 
example  configuration  examined  here,  the  linear  model  should  be  updated  no  less  than  once  per 
4000  cal  of  range.  The  ability  of  lateral  pulsejets  to  change  the  impact  point  steadily  decreases 
with  range  to  the  target. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


^X2^X2 


s~iF  A 
^NA’^NA 


s~iF 

^NPA’^NPA 


/^F  A 
^DD^DD 


/~iF 

^LP^LP 

/~*F  s~\  A 

^MQ^MQ 


IF  IF 

1  XX  9 1  YY 


IA  IA 

1  XX  I1  YY 


rfa’rJy’rfi 


Zero  yaw  axial  force  aerodynamic  coefficient  for  the  forward  and  aft  bodies. 

Yaw  angle  squared  axial  force  aerodynamic  coefficient  for  the  forward  and  aft 
bodies. 

Normal  force  aerodynamic  coefficient  for  the  forward  and  aft  bodies. 

Magnus  force  aerodynamic  coefficient  for  the  forward  and  aft  bodies. 

Roll  moment  aerodynamic  coefficient  due  to  fin  cant  for  the  forward  and  aft 
bodies. 

Roll  damping  moment  aerodynamic  coefficient  for  the  forward  and  aft  bodies. 
Pitch  rate  damping  moment  aerodynamic  coefficient  for  the  forward  and  aft  bodies. 

Viscous  damping  coefficient  for  hydrodynamic  bearing. 

Projectile  characteristic  length. 

Gravity  constant. 

Roll  and  pitch  inertia  of  the  forward  body  projectile  section. 

Roll  and  pitch  inertia  of  the  aft  body  projectile  section. 

Total  projectile  mass. 

Aft  body  mass. 

Forward  body  mass. 

Roll  axis  component  of  the  angular  velocity  vector  of  the  forward  body  expressed 
in  the  fixed  plane  reference  frame. 

Roll  axis  component  of  the  angular  velocity  vector  of  the  aft  body  expressed  in  the 
fixed  plane  reference  frame. 

Components  of  the  angular  velocity  vector  of  both  the  forward  and  aft  bodies 
expressed  in  the  fixed  plane  reference  frame. 

Fixed  plane  components  of  vector  from  composite  center  of  mass  to  forward  body 
mass  center. 
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r  r  r 

ax  5  ay  *  az 


Rfx’RJy’Rfz 


/?/?/? 

’  iVaz 


Rax’Ray’Raz 


Fixed  plane  components  of  vector  from  composite  center  of  mass  to  aft  body  mass 
center. 

Fixed  plane  components  of  vector  from  forward  body  mass  center  to  forward  body 
center  of  pressure. 

Fixed  plane  components  of  vector  from  aft  body  mass  center  to  aft  body  center  of 
pressure. 

Fixed  plane  components  of  vector  from  aft  body  mass  center  to  aft  body  Magnus 
center  of  pressure. 


Rfx ,  Rf ,  Rfz  Fixed  plane  components  of  vector  from  forward  body  mass  center  to  forward  body 
Magnus  center  of  pressure. 

Rix,R1y,Riz  Fixed  plane  components  of  vector  from  composite  body  mass  center  to  the 
impulse  application  point. 


s  Nondimensional  projectile  arc  length. 

x,y,z  Position  vector  components  of  the  composite  center  of  mass  expressed  in  the 

inertial  reference  frame. 

Yj ,  Z;  Pulsejet  force  components  in  projectile  body  axes. 


u,v,w 


V 

6  ,  y/ 


<f> 


F 


f, 


Translation  velocity  components  of  the  composite  center  of  mass  resolved  in  the 
fixed  plane  reference  frame. 

Magnitude  of  mass  center  velocity. 

Euler  pitch,  and  yaw  angles. 

Euler  roll  angle  of  the  forward  body. 

Euler  roll  angle  of  the  aft  body. 
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